The role of domain wall junctions in Carter's pentahedral model 
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The role of domain wall junctions in Carter's pentahedral model is investigated both analyti- 
cally and numerically. We perform, for the first time, field theory simulations of such model with 
various initial conditions. We confirm that there are very specific realizations of Carter's model cor- 
responding to square lattice configurations with X-type junctions which could be stable. However, 
we show that more realistic realizations, consistent with causality constraints, do lead to a scaling 
domain wall network with Y-type junctions. We determine the network properties and discuss the 
corresponding cosmological implications, in particular for dark energy. 



I. INTRODUCTION 

There is now overwhelming observational evidence that 
our Universe is presently undergoing an era of accelerated 
expansion [TJ |2] . In the context of general relativity such 
period can only be explained if the universe is perme- 
ated with an exotic dark energy component violating the 
strong energy condition. The dark energy is often de- 
scribed by a nearly homogeneous scalar field minimally 
coupled to the other matter fields. If the scalar field is 
static then it is equivalent to a cosmological constant but 
the more interesting case is definitely that of a dynamical 
scalar field [3|. 

Nevertheless, the dark energy role is not necessarily 
played by a (nearly) homogeneous field. In fact, it has 
been claimed that a frozen domain wall network could 
naturally explain the observed acceleration of the uni- 
verse @]. However, this possibility has been seriously 
challenged by recent observational results which favor a 
dark energy equation of state parameter, w, very close to 
— 1 (note that w = — 2/3 + v 2 > —2/3 for domain walls, 
where v is the root mean square velocity). Furthermore, 
although it is possible to build (by hand) stable domain 
wall lattices there is strong analytical and numerical evi- 
dence that no such lattices will ever emerge from realistic 
phase transitions [5 , 6j[7||8j[9]. This provided strong sup- 
port for a no-frustration conjecture invalidating domain 
walls as a viable dark energy candidate. 

Still, it has been argued that winding domain wall 
models with X-type junctions could give rise to static 
lattice type configurations thus accounting for at least 
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a fraction of the dark energy density [TOj [TTJ [12j [T3] . 
Carter's pentahedral model JTQI [13] has been constructed 
as an example of a model having an odd number of vac- 
uum configurations giving rise to an even type system 
through the formation of X-type junctions. However, in 
[5j [6] the claim that Carter's pentahedral model would 
form X-type junctions has been challenged and it was ar- 
gued that Y-type junctions would be formed instead. In 
this paper we definitely settle this question. 

Throughout the paper we use units in which c = h = 
m = 1, where the mass scale, m, can be chosen arbitrar- 

iiy. 



II. THE MODEL 

Consider the action 




where <I> and ^ are complex scalar fields, 

C = X + Y -V , (2) 
U($, \I/) is the scalar field potential, 

X = \®*^, (3) 

Y = i**,^, (4) 

and the superscript * stands for the complex conjugate. 
Carter's pentahedral model has a potential given by [TOl 

US] 

V = V ((\t>\ 2 - if + - if) 

+ V e (|$| 2 |*| 2 (cos0 + cos x ) + 2/(1 - e)) , (5) 
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where $ = |$|e^, = |tf|e^, < e < 1, V e = eV , 
= 20 + ip and \ — 2^ ~ 4 > - The potential has five 
minima with V = satisfying cos# = cosx = — 1 and 
|$| 2 = |*| 2 = l/(l-e). 

If e = then Eq. ([5| represents the standard Mexican 
hat potential. Hence, the model allows for cosmic string 
solutions associated with regions where the phase of $ 
and/or \I/ changes by 2n7r, where n is an integer. The 

— 1/2 

string width is roughly S s ~ V Q and consequently its 
energy per unit length is given by 



(6) 



For eCl the static domain wall trajectories connect- 
ing different minima may be calculated to first order in 
e assuming that |3>| = = 1 everywhere. In this case 
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10 ' M ' 10 
with the potential given approximately by 

V(6, X ) = V e (cos 9 + cosx + 2) 

= 2K(cos 2 (<V2) + cos 2 (x/2)) . (8) 

Assuming that |3>| = = 1 everywhere would be 
enough to guarantee that Y-type junctions never occur, 
as long as the energy density remains finite everywhere. 
If the above condition is relaxed then Y-type junctions 
are no longer forbidden and will be associated with cos- 
mic strings (|$| ~ ~ on the string core). 

Consider a planar static domain wall perpendicular to 
the z direction and assume that — 0(z) and cos(x/2) = 
0. The only non-trivial equation of motion is given by 



1 
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or equivalently 



dO 



2F e cos 2 (<9/2), 



cos(0/2) 
which has the solution 



±y/20V e dz, 



tan 



7T + 6> 



(9) 



(10) 



(11) 



for a domain wall located at z = 0. Using Eq. (11) it is 
straightforward to show that 



cot 



=F sinh 



(12) 



where S w = l/y/20V € is the domain wall tickness. In the 
following we shall drop the =p sign. It will be suficient to 
realize that for each solution = 0(z), there will also be 
another solution given by = 6{—z). 




FIG. 1: On the left and right upper panel two possible paths, 
corresponding to Y-type junctions, are illustrated on the sur- 
face of a torus representing the configuration space 
On the left and right middle panel two possible paths, cor- 
responding to stable (left panel) and unstable (right panel) 
X-type junctions, are illustrated. The green line on the right 
middle panel represents a possible decay channel of the un- 
stable X-type junction into two stable Y-type ones. On the 
lower panel the green and red paths (lighter and darker grey 
in black and white) illustrate the trajectories with constant % 
and 0, respectively. 



The energy density, p, associated with the domain wall 



is 



= 4V e cos 2 (6>/2), 



(13) 



where Eq. ([9| was used to obtain the final result. Finally 
the domain wall tension associated with a simple domain 
wall trajectory is given by 



/+oo 
p(z)dz = SV e S w . 
-oo 



(14) 



In the case of a compound wall, in which both and x 
vary along the wall, the corresponding tension is twice 



FIG. 2: Matter-era evolution of a realization of Carter's pen- 
tahedral model with e — 0.2. The simulation starts with ran- 
dom initial conditions with all the five minima having iden- 
tical probability. Note that Y-type junctions are much more 
frequent than X-type ones. From left to right and top to bot- 
tom, the horizon is approximately 1/10, 1/8, 1/6 and 1/4 of 
the box size, respectively. 

that of a simple domain wall 

an = 2(7. (15) 

In Carter's pentahedral model there is a simple domain 
wall trajectory between any of the five minima of the 
potential, with either constant 6 or \. For constant 6 
there is a simple domain wall trajectory between any two 
adjacent minima, with phases ((/>, -0), in the sequence 

1. (7tt/5, -9tt/5) 2. (tt, -tt) 3. (3tt/5, -tt/5) 
4. (7r/5,37r/5) 5. (-7r/5,77r/5) 1. (-3tt/5, 11tt/5) . (16) 

Here and ip vary in five successive steps of ^2tt/5 and 
±4tt/5 respectively thus maintaining = constant (note 
that the phases are defined up to a multiple of 27r) . These 
trajectories are illustrated on the the lower panel of Fig. 1 
by the red path (darker grey in black and white) on the 
surface of a torus with line element 

dl 2 = R\d(t) 2 + (i?i cos + R 2 ) 2 d^ 2 (17) 

where R\ < R 2l representing the configuration space 

If X is a constant then there is a simple domain wall 
trajectory between any two adjacent minima in the se- 
quence 

1. (-137r/5,-97r/5) 4. (-97r/5,-77r/5) 2. (-tt, -tt) 
5. (-tt/5, -3tt/5) 3. (3tt/5, -tt/5) 1. (77r/5,7r/5) .(18) 



FIG. 3: Same as Fig. 1, except that now e = 0.05. 

In this case <p and ip vary by successive steps of ±4tt/5 
and ±2tt/5 respectively thus maintaining x — constant. 
These trajectories are illustrated by the green path 
(lighter grey in black and white) on the lower panel of 
Fig. 1. 

There are Y-type junctions connecting three simple do- 
main walls trajectories. Surrounding a Y-type junction 
there are two domain walls with constant 9 (or x) an d 
another one with constant x ( or #)• A simple example is 
the configuration 

1. (-3tt/5,tt/5) 4. (7r/5,37r/5) 
3. (3tt/5, -tt/5) 1. (7tt/5, tt/5) , (19) 

which corresponds to two domain walls with constant x 
(1 — 4 and 1 — 3) and one with constant (3 — 4). The 
above trajectory is illustrated by the red path on the left 
upper panel of Fig. 1. The overall change in the phase <p 
is equal to 27r. In fact there must always be jump of 27r 
in either <p or ip around a Y-type junction. 

Another example corresponding to a Y-type junction 
is the configuration 

1. (-3tt/5,7t/5) 2. (-tt,7t) 
5. (-tt/5, 7tt/5) 1. (-3tt/5, 11tt/5) , (20) 

where two domain walls with constant (1 — 2 and 1 — 5) 
and one with constant x (2 — 5) meet. In this case it is the 
overall change in ip that is equal to 2tt. This trajectory 
is illustrated by the red path on the right upper panel of 
Fig. 1. 

What about X-type junctions ? Is there a trajectory in 
which <j) and ip are continuous around a X-type junction 
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? The answer is yes. For example, both <fi and ip can be 
made continuous around the X-type junction described 
by the following configuration 



1. (77r/5,7r/5) 2. (7r,7r) 4. (7r/5,37r/5) 
3. (3tt/5, -tt/5) 1. (7?r/5, tt/5) . 



(21) 



This trajectory is illustrated by the red path on the 
left middle panel of Fig. 1. As correctly pointed out 
in [TQl [13], Carter's penthedral model allows for square 
domain wall lattice solutions that are stable, if e is suffi- 
ciently small. However, as we will show in the following 
section, such lattices are never generated from realistic 
initial conditions. 

Around a X-type junction where three walls with con- 
stant (or x) meet one wall with constant \ (or both 
(j) and ip must change by a factor of 2tt. Consider the 
following example which is illustrated by the red path on 
the right middle panel of Fig. 1 

1. (-37r/5,-97r/5) 2. (-7r,-7r) 5. (-tt/5, -3tt/5) 

3. (371-/5,-71-/5) 1. (7tt/5,7t/5). (22) 

In this case the energy of the junction associated with the 
presence of a string is greater, by a factor of 2, compared 
to Y-type junctions. Hence, the string does nothing for 
the stability of the junction. Such X-type junction would 
be unstable and decay into a pair of Y-type junctions, 
even if e is small (the green line represents a possible 
decay channel). This is the reason why, in the context 
of Carter's pentahedral model, Y-type junctions are pre- 
ferred, with exception of very specific realizations. 



III. SIMULATIONS 

In order to test our analytical expectations we will now 
present the results of a few 256 2 simulations in two spatial 
dimensions. Although these simulations are relatively 
small in size and dynamical range, they are more than 
enough to support our analysis. In all the simulations we 
use the PRS algorithm [14] modifying the domain wall 
thickness in order to ensure a fixed comoving resolution. 
More details about the numerical code can be found in 
[9] and references therein. 

Fig. 2 shows four snapshots of a matter era simulation 
of a realization of Carter's pentahedral model (e = 0.2) 
with random initial conditions. At each grid point, one 
of the minima was randomly assigned, all the minima 
having equal probability. The cosmic time, t, is increas- 
ing from left to right and top to bottom (the horizon is 
approximately 1/10, 1/8, 1/6 and 1/4 of the box size re- 
spectively). The simulations show that Y-type junctions 
are much more frequent than X-type ones. This is not 
surprising since the probability that the combination of 
two Y-type junctions will give rise to one stable X-type 
can be easily calculated and is equal to 2/9, assuming 
that the corresponding minima are randomly chosen with 



FIG. 4: The configuration space distribution for the last time 
step of the simulation in Fig. 2. On the left panel the x and 
y axis represent the phases <j> and -0, respectively. On the 
right panel the x and y axis represent Re(<E>) = |<I>| cos0 and 
Im(<3>) = |<3>| sin0. 





FIG. 5: The same as in Fig. |3]for e = 0.05. 



equal probability, subject to the constraint the same min- 
ima cannot be assigned to both sides of a domain wall. 
On the other hand, the probability that the collapse of 
a square domain with Y-type junctions at the vertices 
will give rise to a stable X-type junction is equal to 1/21, 
again assuming a random configuration. Furthermore, 
this does not take into consideration that stable X-type 
junctions may break into two Y-type junctions if enough 
energy is available. However, some rare stable X-type 
junctions can still be identified in the simulations. 

Fig. 3 is similar to Fig 2 except that now e = 0.05. As a 
consequence, the energy density inside the domain walls 
is reduced by a factor of 4 while their thickness increases 
by a factor of 2. On the other hand, the strings remain 
roughly the same. Of course, in the e — > limit the dy- 
namics would be completely dominated by the strings. 
However, in this limit the thickness of the domain walls 
becomes very large (5 W oc e -1 / 2 ) and the domain wall 
network would no longer be well defined. In any case, 
this would not help domain walls as a possible dark en- 
ergy candidate since, in that case, the contribution of the 
junctions to the energy density would be the dominant 
one, thus leading to an equation of state parameter signif- 
icantly greater than —2/3. Moreover, the strings have a 
small impact on the overall dynamics as long as the aver- 
age domain wall energy density dominates over that asso- 
ciated with the junctions. This happens for aL ^> /i ~ 1 
or equivalently S w /L <C 1, where L is the characteristic 
scale of the network. Such condition is always verified as 
long as the thickness of the domain walls is much smaller 
their typical curvature scale. In fact, the string energy 
per unit length (/i ~ 1) of a stable Y-type junction is of 
the same order as the energy per unit length of a stable 
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FIG. 6: Evolution of a periodic square lattice of Carter's pen- 
tahedral model. The initial configuration of minima was cho- 
sen to allow for both stable and unstable X-type junctions. 

X-type junction (~ V e S^ ~ 1), which means that X-type 
junctions are configurations of delicate equilibrium, sus- 
ceptible to decay in the presence of relatively small per- 
turbations. Hence, even the rare stable X-type junctions 
which appear in the simulations would probably not be 
there if the domain wall thickness had not to be artifi- 
cially enlarged in order to ensure that the domain walls 
were resolved by the numerical code. 

Fig. 4 shows the configuration space distribution for 
the last time step of the simulation in Fig. 2. On the left 
panel x and y axis represent <p and iJj : respectively. The 
left panel shows that only simple domain wall trajecto- 
ries with constant — 2</> + ip (or x — 2^ — 4>) appear 
in the simulations. On the right panel, the x and y axis 
represent Re(<I>) = |<£|cos(/> and Im(<3>) = |<3>|sin0, re- 
spectively. The five different minima corresponding to a 
constant value of |3>| = yl/ 0- ~ e ) — 1-H? as we ^ as 
the corresponding domain wall trajectories in <l>, can be 
easily identified. Fig. 5 is similar to Fig 4 except that 
now e = 0.05 so that the minima on the left correspond 
to | $ | ~ 1.01. As a result, the domain wall trajectories 
appear as nearly circular orbits on the left panel of Fig. 5. 

Fig. 6 shows the evolution of a hand-made periodic 
square lattice realization of Carter's pentahedral model 
with e = 0.2. The initial configuration of minima was 



chosen to allow for X-type junctions corresponding to a 
continuous (j) and and X-type junctions around which 
both (j) and ip change by a factor of 2tt. As expected, 
the simulations show that the former are stable while 
the later are unstable and decay into two stable Y-type 
junctions. 

It is also possible to choose the initial conditions in 
a way that a square lattice with only X-type junctions 
is formed and we have verified that such a configuration 
is stable, as claimed by Carter [TUl [13]. However, one 
should bear in mind that it corresponds to a very specific 
set of initial conditions which would violate causality, if 
they were to extend over scales larger than the particle 
horizon. 



IV. CONCLUSIONS 

In this paper we confirmed that there are very special 
realizations of Carter's pentahedral model, corresponding 
to square lattice configurations with X-type junctions, 
which could be stable. However, we have shown that 
more realistic realizations of Carter's pentahedral model, 
such as those with random initial conditions, give rise to 
a network with Y-type junctions. This leads to a domain 
wall network whose properties are virtually indistinguish- 
able from those of a specific realization of the ideal class 
of models with 4 real scalar fields (and 5 minima), with 
similar initial conditions. The ideal class of models has 
been studied in detail in [7J [9] where a compelling evi- 
dence for a gradual approach to scaling, with L oc t, was 
found both in the radiation and matter eras. As a result, 
and in spite of its very interesting topological properties, 
Carter's pentahedral model does not naturally lead to a 
frustrated network with v ~ and L « t, a necessary 
condition for domain walls to provide a contribution to 
the dark energy budget. There are other models which 
allow for X-type junctions (see for example, [15j[T6]) but 
they also do not lead to a frozen network, starting from 
random initial conditions |7J |9] . 
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